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ABSTRACT 

An advanced 3-D Computational Fluid Dynamics (CFD) model was developed to analyze the 
flow interaction between a gas turbine combustor and an integral bleed plenum. In this model, the elliptic 
governing equations of continuity, momentum and the k-e turbulence model were solved on a 
boundary-fitted, curvilinear, orthogonal grid system. The model was first validated against test data from 
public literature and then applied to a gas turbine combustor with integral bleed. The model predictions 
agreed well with data from combustor rig testing. The model predictions also indicated strong flow 
interaction between the combustor and the integral bleed. Integral bleed flow distribution was found to 
have a great effect on the pressure distribution around the gas turbine combustor. 

INTRODUCTION 

The onboard gas turbine auxiliary power unit (APU) generally is designed to deliver bleed air 
either from a basic APU or from a separate compressor. The bleed air supplied by the APU normally is 
used for starting the main engines and operating the air-cycle air conditioning and pressurization system. 
On the direct-bleed APU (also called integral-bleed APU), bleed air is extracted from a point between the 
compressor and combustor, as shown in Figure 1, and is routed to the airplane pneumatic system via a 
bleed valve attached to a port on the turbine plenum. 

The bleed air extraction from the turbine plenum has long been recognized among APU 
combustor designers to have a significant effect on combustor performance. However, knowledge in this 
field is limited and has been acquired only as a result of rig or engine testings. No theoretical modeling 
has been conducted and reported in the literature. The purpose of this study is to do a comprehensive 
theoretical investigation into the flow interaction between a gas turbine combustor and an integral bleed 
plenum. Findings from this study should prove useful in the design and development of a gas turbine 
combustor with integral bleed. 

The flow interaction between a gas turbine combustor and an integral bleed is fully 
three-dimensional and extremely complicated. To do 3-D Computational Fluid Dynamics (CFD) 
modeling of this complex flow interaction using a cylindrical grid system requires a huge number of grid 
nodes and complex treatment of wall boundary conditions. Recently, a large variety of numerical grid 
generation methods [1] have been developed to simplify wall boundary condition treatment for any 
complex flow geometry. These methods can generate either orthogonal or non-orthogonal 
boundary-fitted grid systems for any complex flow geometry. In this study, a boundary-fitted, 
curvilinear, orthogonal grid system is used for the two-dimensional (x,y) plane and the angular (z) 
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Figure 1. A Typical Gas Turbine Combustor with Integral Bleed. 


dimension is considered as a body of revolution. The choice of a curvilinear orthogonal grid system in 
the present flow modeling is motivated by easier wall boundary condition treatment, inclusion of most of 
the grid nodes in the computational domain, and no extra terms of cross-derivative appearing in the 
transformed governing equations. 

ORTHOGONAL GRID GENERATION 

As mentioned above, the flow geometry of a gas turbine combustor with integral bleed is quite 
complex for 3-D CFD modeling. A boundary-fitted, curvilinear, orthogonal grid system is therefore used 
here. The orthogonal grid system used in this study was obtained numerically based on the following 2-D 
orthogonal grid generation model: 
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Governing Equations: 
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where (x,y) are the Cartesian coordinates of the grid points in the physical plane, as shown in Figure 2, 
and (S,n) are the coordinates of the corresponding grid points in the transformed plane. The control 
functions P and Q in equations (1) and (2) are used to concentrate grid lines as desired. The function F 
denotes user-specified boundary curves. In order to facilitate numerical solution for the above equations 
(1) to (6), they were transformed into the following equations: 

Governing Equations: 

ax^ - ♦ Y^„ * + x n Q) (la) 
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Figure "L Coordinate Transformation. 
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where 
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Boundary Conditions: 

SCI: + y<y„ » 0 ; G t U,y) =0 (3a) 
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These transformed equations are solved numerically using a finite difference technique. This grid 
generation modei only generates a boundary- fitted, curvilinear, orthogonal grid system in two 
dimensions. The present flow calculation requires a 3-D orthogonal grid system, obtained by rotating this 
2-D orthogonal grid system about an axis in the same plane. Using this approach, a typical 3-D grid 
system for the combustor with integral bleed flow calculation is shown in Figure 3. 

MATHEMATICAL MODEL 

The transport equations for a 3-D incompressible turbulent flow can be written as 

- r§> . - r|&> ♦■£<?*- rjf) - * 

where represents a general variable, p is density, U, V, and W are velocities, r is effective turbulent 
diffusivity, and 5* is the source term for variable $ . 

The above transport equations are given in Cartesian coordinates (x,y,z). Prior to their numerical 
solution, these equations are transformed into general orthogonal coordinates ( 5 , n / C )• The resulting 
transport equations are expressed as follows: 

- ■£■§£>) - 

where h represents scalar coefficients. 

In order to solve the above partial differential equations governing the flow of fluid through a 
combustor annulus with integral bleed, boundary conditions must be specified at the appropriate 
locations. In the present study, the boundary conditions at the annulus inlet and the bleed duct exit are 
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provided from the engine cycle analysis, while the airflow through the various orifices in the combustor, 
as shown in Table 1, are calculated using a I-D annulus flow model. 
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Figure 3. A Typical Orthogonal Grid System for the Combustor 
with Integral Bleed Flow Calculation. 
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TABLE 1. AIRFLOW DISTRIBUTIONS AROUND THE COMBUSTOR 




on 

H2 

Primary Holes 


4.28 

9.90 

Dilution Holes 


5.59 

5.69 

Film Cooling 


10.90 

7.15 

Dome Cooling 


11.21 


Bleed 


42.01 


Shroud Air 


1.10 


Leakage 


J.28 

ois> 

Total 


76.37% 

23.63% 


The solution of the curvilinear transport equations is accomplished using a finite volume, 
structured grid, sequential solution, numerical method. The differential equations are integrated over 
discrete volumes in the domain of interest using assu mptio ns of linear and stepped profiles to obtain their 
numerical counterparts. Interaction between the convective and diffusive terms is handled numerically 
using Upwind Hybrid Differencing [2], while the coupling between the continuity and momentum 
equations is treated by the SIMPLER algorithm [3] and the turbulence closure is accomplished by the 
standard k-e model. At each iteration step, the set of linear equations is solved using the Whole Field 
Solver (WFS) routine of Przekwas [4]. The WFS has proven to be very efficient as it is formulated for 
structured grids and cyclic boundary conditions which are used in this analysis. 

The use of orthogonal coordinate systems adds effort to the grid generation process, as numerical 
equations must be solved to determine the position of the coordinates in the physical space. It is also 
frequently the case that the shape of the domain is such that a completely orthogonal grid system cannot 
be generated and a certain amount of deviation from orthogonality must be tolerated. However, the 
benefits of this extra work are realized when the flow equations are solved. The terms associated with 
grid non-orthogonality are usually treated in an explicit manner and in situations of skewed control 
volumes, can be the same order of magnitude as the normal implicit convection and diffusion terms. As a 
result, convergence can be hindered. Orthogonal grids create a more implicit nature to the numerical 
equations and have inherently better convergence rates. 

The combination of the SIMPLER algorithm and the orthogonal grid structure has proven to be 
extremely robust. In all but a very few situations, relaxation factors of 0.5 are used on all equations 
(except for pressure and the velocity correction which require 1 .0 relaxation factors). Converged 
solutions using up to 400,000 cells have been obtained in the order of 750 iterations. At convergence, the 
sum of the absolute continuity error of all the control volumes, normalized by the total system flow rate, 
is less than 0. 1 percent. 


NUMERICAL RESULTS 

The predictions of the CFD model used in this study have been compared to numerous sets of 
experimental data found in the literature. In order to illustrate the predictive capability of the code, two 
such cases have been included below, both of which are 2-D flows. For these situations, the 3-D code 
was run with three circumferential planes. 

The first case is from the experiment of Roback and Johnson [5] and consists of two coaxial 
flowing water jets that discharge into a confined region (illustrated in Figure 4). The inner jets are 
swirling while the outer are not. Velocity data was taken at several axial stations downstream of the jet 
discharge. Figure 5 shows the comparison of the predicted and measured axial, radial, and tangential 
velocity components respectively. The comparison is as good as should be expected considering the well 
known deficiencies of the k-e model for swirling flow and the uncertainty in the boundary conditions of 
the data. 
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Figure 4. Axisymmetric Swirling Flow Geometry. 


The second case described is that of Sovran and Klomp [6] and consists of a 6-degree wall 
annular diffuser with an area ratio of 1.955. The computed static pressure contours (Pascals) are shown in 
Figure 6. The experimental data for this case consist of the static pressure recovery coefficient, Cp. For 
the configuration analyzed, the measured Cp was 0.6 which should be compared to the value of 0.643 (7 
percent high) calculated from the CFD model output, indicating resonable accuracy in the pressure 
calculation methods used in the code. 

After validation of the 3-D CFD model using experimental data found in the literature, the present 
model was then used to simulate complex external flow of a gas turbine combustor with integral bleed. 
The computed results using a grid system of 91x21x97 nodes, as shown in Figure 3, are discussed in the 
following: 

Figure 7 shows the predicted U-V velocity field at K-planes 1, 38, and 59. The predicted results 
illustrate a strong recirculation zone formed near the annulus inlet. The formation of this recirculation 
zone is due to the presence of a backward-facing step in the flow path. This recirculation zone has caused 
pressure loss near the combustor O.D. dilution holes. As a result, the dilution jet penetration will be 
reduced because of less pressure drop across dilution holes. Based on gas turbine combustor design 
experiences, inadequate dilution jet penetration usually cannot break up the hot gas regions efficiently and 
will result in high pattern factor at the combustor exit. 

The predicted results, as shown in Figure 7, also illustrate non-uniform bleed air extraction from 
the turbine plenum. The bleed air extraction at K-plane 59, which is closer to bleed port (K-planes 
61-65), is higher than that at K-plane 38. The non-uniform bleed air extraction indicates that the present 
integral bleed design has not been optimized and has caused some flow separations in the annulus. It is 
also interesting to note that bleed air enters the bleed duct like an impinged jet and creates a vortex pair in 
the duct. 

Figure 8 shows the predicted U-W velocity field on the unwrapped surface of combustor liner. 
The predicted results indicate that a forward stagnation point, where the air is brought to rest with an 
accompanying rise in pressure, was formed upstream of a fuel atomizer or an ignitor. The predicted 
results also indicate that a strong wake was formed downstream of a fuel atomizer or an ignitor. 

The presence of eight fuel atomizers and one ignitor has a profound effect on pressure 
distributions around the combustor liner. This can be seen in Figure 9, where the predicted static pressure 
distributions around the combustor liner is plotted using color graphics. In this figure, each fuel atomizer 
was found to be associated with a high pressure region upstream and a low pressure region downstream. 
The same pattern was observed for the ignitor. Figure 9 also shows non-uniform bleed extraction has a 
great effect on pressure distribution around the combustor liner. This is illustrated by results of 
non-uniform pressure distribution on the dome. Based on gas turbine combustor design experience, 
non-uniform pressure distribution around the combustor liner usually causes non-uniform airflow 
distribution and will result in hot and cold spots at the combustor exit. 

Figure 10 shows the predicted static pressure distributions around the turbine plenum. The 
predicted results further illustrate that the presence of eight fuel atomizers and one ignitor has a 
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Figure 6. Pressure Contours for Annular Diffuser. 
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Figure 7. Predicted U-V Velocity Field around Combustor. 
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Figure 9. Predicted Pressure Distributions around Combustor. 









significant effect on pressure distributions around the turbine plenum. The predicted results also show 
non-uniform pressure distributions on the bleed duct wall. The predicted results of non-uniform pressure 
distributions on the bleed duct wall (at 1=43 and J=20) were compared to combustor rig test data in Figure 
1 1. This comparison indicates model predictions agreed well with test data from combustor rig. 

Having obtained good predictions of non-uniform pressure distribution on the bleed duct wall, it 
is then possible to consider another 3-D case run in order to study the effect of intergral bleed extraction 
on the pressure distributions around the combustor liner. Figures 12 and 13 show the predicted pressure 
distributions at various hole locations for a gas turbine combustor with and without an integral bleed. 
Comparison of Figures 12 and 13 confirms that integral bleed flow distribution has a great effect on the 
pressure distributions around the gas turbine combustor. 
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Figure 11. Comparison of Predicted and Measured Pressure 
Distributions on the Bleed Duct Wall. 


PUSCtDINS Puns ULUHK NOT KK.MED 


377 


STATIC PRESSURE, PSI STATIC PRESSURE, PSI STATIC PRESSURE. PSI 


74.00 -i 


OD - DO (1-18) 


74.00 -i 


DOME -ID (I -51) 


73.60 A 



72.00 A 1 1 1 1 

0 90 180 270 360 

ANGLE FROM TDC, DEGREES 


£ 73.60 - 
Ui 

1 73-20 - 

& 

ol 7Z80 - 
O 

5 

K7Z40- 

72.00 -I 1 1 1 n 

0 90 180 270 360 

ANGLE FROM TDC. DEGREES 





74.00 
73.60 A 


73.20 -i 


7Z80 -I 


72.40 A 


72.00 


DOME-OD (I -39) 




1 1 ! 1 

0 90 180 270 360 

ANGLE FROM TDC, DEGREES 


74.00 n 


£ 73.60 


UJ 

<r 

</3 

w 

Ui 

QC 

a. 

o 


73.20 H 


72.80 H 


S 7Z4 ° H 


72.00 


ID - DO (I » 73) 


Mvvvwv wwvwWVWMl 


1 1 1 r 

0 90 180 270 360 

ANGLE FROM TDC, DEGREES 


GC11475-t2 


Figure 12. Predicted Pressure Distributions 
at Various Combustor Hole Locations. 
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Figure 13. Predicted Pressure Distribution at Various 
Combustor Hole Locations (without Integral Bleed). 
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CONCLUSIONS 


An advanced 3-D CFD model based on a boundary-fitted, curvilinear, orthogonal grid system has 
been developed. The model is capable of predicting the complex external flow field of a gas turbine 
combustor with integral bleed. Several important findings from this study of Combustor-Integral Bleed 

flow interaction are summarized below: .... 

(1) Integral bleed flow distribution has a great effect on the pressure distribution around the gas 

turbine combustor. . , 

(2) The presence of atomizers and ignitors has a profound effect on pressure distribution aroun 

the combustor O.D. liner. Each atomizer or ignitor was found to be associated with a high pressure region 
upstream and a low pressure region downstream. 

(3) The non-uniform pressure distribution usually causes non-uniform airflow distribution around 

the combustor liner and will result in hot and cold spots at the combustor exit. 

(4) The recirculation zone formed near the annulus inlet has caused pressure loss near the 
combustor O.D. dilution holes. As a result, the dilution jet penetration is reduced and will result in high 
pattern factor at the combustor exit. 
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